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A METHOD FOR DETERMINING AERODYNAMIC DAMPING 
OF THE STRUCTURAL BENDING MODES WITH APPLICATION 
TO THE APOLLO/SATURN I LAUNCH VEHICLE 

SUMMARY 


Calculation of the aerodynamic damping of the Apollo/Saturn I 
Launch Vehicle has been made for the first two structural bending modes. 
A comparison of the damping derivative rith wind tunnel data indicates 
that this method is valid for analytic determination of the aerodynamic 
damping for a vehicle of this type. 


INTRODUCTION 

Air forces due to the n^ modal velocity are generally negligible 
compared to structural design airloads but provide a beneficial source 

of damping on the n 1 mode itself. In the ease of a finned body such 
as the Apollo/ Saturn I the aerodynamic damping can be larger than the 
structural damping and of considerable importance. 

The modal velocity induces small changes in angle of attack along 
the vehicle. The distribution of lift due to the changes in angle of 
attack may be determined by slender body theory (with some modification 
for the command module). From this lift distribution, the generalized 
force is determined for the modal differential equation. The ratio of 
aerodynamic damping to critical damping as a function of Mach number is 
found from this equation. 

This paper presents a technique for obtaining the aerodynamic 
damping from the modal velocity. Application of this technique to the 
Apollo/ Saturn I configuration consists of calculations made for a 
specific time in the flight trajectory and extension of these calcula- 
tions to other times of interest. 
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y n (x,t) 


Normal force coefficient per unit angle of attack, 
1/deg 

Generalized force associated -with mode n. I'd 


Vehicle length, in. 
2 


Mass. 


lb-sec 


ft 

Generalized mass 

Normal force, lb 
Mode number 

2 

Free-stream dynamic pressure lb/ft 

2 

Reference area of vehicle as a function of x, in. 

o 

Time, sec 
Velocity, ft /sec 

Distance along vehicle (zero at the nozzle end), in. 
Deflection in mode n 


a 

P 

P 
(0 

Subscripts 

A Aerodynamic 

cr Critical 


Angle of attack, radians 
Damping coefficient 

3 

Stream density, slug/ft 
Modal frequency, rad/ sec 


s 


Structural 


Modal number 


' 1.0 


ANALYSIS 
Method 

Sketch (a) presents the 
Apollo/Saturn I vehicle as 
deflected in the first free- 
free bending mode. Displace- 
ments and velocities in each 
structural mode induce aero- 
dynamic forces on the vehicle 
due to local angle of attack 
along the vehicle. For modal 
displacement _ dy ^ 

a » S . 

dx 

The aerodynamic force distribution due to 



<ty n (x) 

dx 


Sketch (a) 


for slender missiles 


is sm al l and does not affect damping; therefore, it is neglected for 
this vehicle. 


The change in local angle of attack 


a 

x 


due to modal velocity 


produces normal force loads that induce aerodynamic damping. Although 
these loads are small for vehicles such as the Apollo/Saturn I, they 
can produce a significant amount of damping. This damping is dependent 
on the vonicle configuration, modal shape and modal frequency* A 
drawing of the Apollo/Satum I launch vehicle is given in figure 1. 

The local angle of attack a due to modal velocity can be 
approximated, provided a is small, by: 




n 


(i) 



The modal displacement when normalized can be written as: 


y n (x,t) 


y n (x)sin a> fl t 


( 2 ) 


k 


so that 

y n (x,t) = a> n y n (x) cos a> n t 


For t = 0 or maximum value of y^ 


*„<*> ‘ “n y n (x> 


( 5 ) 


Thus, substituting (5) in (l) 



- 0 ) y (x) 
n ijr 

V 




The distribution of lift for small changes in a caused by y^(x) may 


be determined by use of slender body theory for bodies of revolution 

(ref. 1) 

N * 2a qS 
n 

(5) 

so that 

^(a n S) 


or 

« . - 24 "n a , 3, 

dx V dx * y n S ' 

(6) 

or 

• 



S * - " “n £ < y n S > 

(6a) 


This equation shows that the lift due to modal velocity is a function 
of the forward velocity of the vehicle instead of the square of the 
velocity, as expected. 


Lift Calculations 

In order to calculate this normal force along the vehicle, parameters 
from the 8A*5 trajectory, as given in reference 2, were used with modal 
shapes and frequencies calculated in a digital program for a particular 
time* She time selected was for m a xim um dynamic pressure, 67 seconds 
after lift off* The dynamic pressure was multiplied hy 1.2 i# orderja, . 
take into f$eouot possible deviations from the nominal traj$c$<^. '* 
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first two bending modes normalized at the nose are shown plotted in 
figures 2 and 5. The important data associated with this case for 
t = 67 seconds are: 


Mach No. * 1.55 

q = 864 lb/ft 2 
V « 1,400 ft/sec 
=11,32 rad/sec 
ajg = 23.05 rad/sec 

The lift distribution is shown plotted in figures 4 and 5* The data 
points shown on these figures are from wind tunnel tests on the Command 
Module and were used for the lift distribution. The generalized force 
due to modal velocity is 



(7) 


In order to determine the damping due to unit modal velocity, the 
generalized force must be divided by the modal frequency giving 
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Calculations for the first two modes were made to determine 
the time t * 67 seconds. 


( 8 ) 



for 


The first two modes for t » 67 seconds given in figures 2 and 5 were 
assumed to be constant over the Mach range of interest, M » 0.7 to 

M * 1 . 55 . Using this assumption was calculated for these Mach 

c ■*» 

numbers by using the ratio of ^ along the flight trajectory to the 
value at M • 1.55* The generalized forces due to unit modal frequency, 
excluding fins are: 


<1© 


= - 1.0869 


and 



The lift for the fins due to local angle of attack change was obtained 
.separately using the C for M = I.'Xj obtained from wind tunnel tot l 

X * a n 

(ref. j ) . This value was determined by taking the difference in 0^ 

on the vehicle with fins and C on the vehicle without fins at M - 1. 

n 

For M = 1.55 

= 0.057 r~~ (Based on 3 * /oO ft 2 ) 
n a *" deg 


Center of pressure of fin = 175 inches from station 0 


a = 
n 


-n y n< 175) 


Oj, = — = -0.00000064 (1st mode) 
a 2 « " -■ ^§^ 5 ^ = 0.000192 ( 2 nd mode) 


Dividing through by co^, the concentrated lift forces on the fin become 



(.056H57.3)(»0.00006qC4)(864)(500) 

11.52 


and 


- -5.57 


lb 

rad/ sec 


"a . (. 0561(37. 3)(O.OOQ193)(86k)f5CO) 

<% 35.03 


8.28 — 

• rad/sec 


« 



Lift due to the fins is given as concentrated force at the fir. center 
of pressure in figures 4 and 5* Generalized forces due to unit modal 
frequency on the fin: 


4 


Thus 


and 



y l 

~ = -?o7 (0.093) 

j. 

= -O.5O3 

— = 3.23 (-0.14) 

= -1.1592 


Total generalized forces due to unit modal frequency: 


and 




%( along vehicle) 


+ 



o> 


n(fins) 



-1.241;, -0.503 
-1.7V-5 



-1.0809 -1.1^92 


2.244 
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A form of the generalized equation of motion for elastic modes is: 


2 * n , 

My = -o> M y T y 

n n n n 'n cd n 

n 


U 

where M = generalized mass = / M(x)y 2 (x)dx 

n yj 


and M(x) = distributed mass 

Thus, the generalized equations become: 

18.708 y ± + 1.746 y x + ( 11 . 52) 2 (18.708^ = 0 (1^ = 18.708) 

and p 

22.11 y 2 + 2.244 + (23.04)^ (22.1l)y 2 *0 = 22.11) 

The ratio of aerodynamic damping to critical damping is shown in the 
following differential equation (ref. 4) : 

F. F. 

y 2 

My- — — y + a) M y = 0 Let P — * coefficient of 

n n <0 n n n n o> . . . „ 

n n damping term 


In operator form 

2 fi 2 

D + „ D + o> *0 
w n 

n 




whose solution is: 


/^. 4o> 2 
2 n 




+ K„e 


I 



or 
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The behavior of the damped system depends upon whether the radical 
of this equation has a value which is real, imaginary, or zero. For a 
critical damping coefficient P the value of the radical must be zero. 


or 


0 * 2 a> M 

cr n n 


The amount of damping in the system is then determined by the ra_:io 


0/Pc 


J_ = _ 


2o> M 
n n 



2(o M 
n n 


For the first mode, this becomes 



The second mode gives 
B 2,244 

P^ “ 2(23.05) (22.109) = °‘° 021 

These values are jotted in figure 6 for M = 1.55> along with the 
values calculated for the Mach numbers considered. 


The values for the other Mach numbers were obtained by multiplying 
the generalized force along the vehicle value for M = 1.55 by the 
ratio 


qA Mach No. « 0.7, 0.8 > . . . . 1.5 
q/V Mach Ho/« 1.5 


The generalized force values for the fins were obtained in the same 
manner ; however, the variation of C n of the fins with Mach number 
was also taken into account. a 
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RESULTS 


The ratio of aerodynamic damping to critical damping is presented 
in figure 6. Also in this figure are the results from the Langley 
SD-1 shake test on the Apollo/Saturn I configuration (unpublished data) . 
The following results may be noted: 

1. The structural damping (as indicated by the dashed line) is 
only a small percentage of the total damping for both the calculated 
values and the tests results for the model investigated. 

2. An increase in the calculated aerodynamic damping at Mach 

Q 

number 1.0 is due primarily to an increase of the n of the fins in 

this area. This increase in damping does not appear in the experimental 
results. The overall difference in magnitude of the calculated and 
experimental results may be attributed in part to the difference In 
trajectory parameters and modal shapes. 

3. The result of a study on separated flow effects on the 
Saturn I is given in reference 5. This study indicated that the 
separated flow along the vehicle has a large effect on the steady and 
unsteady loads. Calculations using these loads give a larger percent 
of aerodynamic damping them those plotted in figure 6, particularly in 
the second mode. Part of this increased damping is due to the large 
positive load on the Command Module induced from the wake of the disk 
on the escape rocket. Since the disk has been removed from the tower 
configuration there would be a corresponding decrease in the aerodynamic 
damping as determined in reference 5* The remaining difference cannot 
be fully accounted for without a review of the calculations used to 
determine the damping in reference 5* 

CONCLUDING REMARKS 

The method used appears to be a valid approach for the determination 
of aerodynamic damping on the Saturn I. Damping in the first mode due 
to the aerodynamic effect is of a much larger magnitude than the struc- 
tural damping, with the largest amount being at Mach number 1.0. The 
second mode damping calculations show a smaller ratio of aerodynamic 
damping to structural damping; however, this amount would still consti- 
tute appreciable damping in the second mode. 
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Figure 2 .- First mode deflection (t * 67 sec), col = 11. 32 rad/sec. 





67 sec). <cu » 25*05 rad /t 





Figure 4.- Lift due to first mode velocity. 
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Figure 5»- Lift due to second mode velocity 






Variation of damping with Mach number 
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Figure 6.- Concluded. 





